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Abstract. We present mid-IR imaging observations of a high Galactic latitude cirrus obtained with the ISO camera 
ISOCAM at 6" angular resolution. The observations were done with two filters LW2 (5-8.5 /im ) and LW3 (12-18 
fim ) that measure respectively the aromatic carbon bands and the underlying continuum emission from small 
dust particles. Three 0.05 square degree images sample atomic and molecular sections in the Ursa Major cirrus. 
These images are compared with Hi , CO and IRAS observations. In such a cloud transparent to stellar light 
(Av < 0.5) the mid-infrared to 100 fira and the mid-IR emissivity per hydrogen are related to the abundance and 
the optical properties of small dust particles independently of any modelling of the penetration of the radiation. 
Within the atomic section of the cloud, the comparison between ISOCAM images and 21 cm interferometric data 
highlights an enhancement of the mid-IR emitters abundance by a factor ~ 5 in an Hi filament characterized 
by a large transverse velocity gradient suggestive of rotation. Furthermore, a drop in the abundance of the same 
mid-IR emitters is observed at the interface between the atomic and molecular cirrus sections. We propose that 
these abundance variations of the mid-IR emitters are related to the production of small dust particles by grain 
shattering in energetic grain-grain collisions generated by turbulent motions within the cirrus and inversely by 
their disappearance due to coagulation on large grains. At the H1-H2 interface we also observe a change in the 
J„(LW2)//„(LW3) ratio by a factor 2. This color variation indicates that the amplitude of the continuum near 15 
yitm, relative to the aromatic bands, rises inside the molecular region. It could result from a modification of the 
dust size distribution or of the intrinsic optical properties of the small dust particles. 

Key words, interstellar medium; infrared cirrus; dust; turbulence. 



1. Introduction 

We know that interstellar dust plays an important role in 
the physics of the interstellar medium. For instance, dust 



grains heat th e gas via the photoelectric effect (Bakes & 
Tielens, 1994 ) an d provide the surfaces nec essary to the 
formation of H2 ( Duley fc Williams, 1993] ). These pro- 
cesses are sensitive to the size distribution, the struc- 
ture and composition of dust grains and, especially, to 
the abundance of the smallest particles. The IRAS images 
of nearby molecular complexes and bright cirrus clouds 
(Ay > 1) have been extensively used to compare the spa- 
tial distribution of the mid-IR emission from the smallest 
stochastically heated grains and that of the far-IR emis- 
sion from larger grains in thermal equilibrium. The color 
ratio between these two emissions has been related to 
the abundance ratio between small and large dust grains 
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which was inferred to widely vary from cloud to cloud and 
within clouds (Boulanger et al., 199C). The gain in bright- 
ness sensitivity and angular resolution provided since then 
by dedicated observations carried out with ISOCAM, the 
camera on board the Infrared Space Observatory (ISO), 
now allows us to carry on these IRAS investigations and 
in particular to extend them to the range of column den- 
sities (Nh ~ 5 x 10 20 cm -2 ) where the transition from 
atomic to molecular gas takes place. 

ISOCAM has provided low resolution spectra (A/ AA ~ 
50) of the mid-IR interstellar emission over a wide range 
of environments including the diffuse Galactic emission 
and cirrus clouds. The emission in the 3-15 /im range is 
invariably associated with emission bands considered to 
be characteristic of aromatic hydrocarbon particles. To 
be heated to temperatures of a few 100 K by single stel- 
lar photons these particles have to be smaller than a few 
thousands carbon atoms. The exact nature, structure and 
origin of the emitting particles are still a matter of debate. 
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Fig. 1. The Ursa Major cirrus observed by IRAS at 
100 fim (color map) with the CO integrated emission of 
Pound fc Goodman (1997 ) in contours (first level = 5 K 



km s I and slep is alsu 5 K kin s — ). The pusiliuu of Lhe 



three IBOGAM fields ubserved are alsu shuwu. 

The abundance variations revealed by IRAS images show 
that small dust grains might not be related to the parti- 
cles released by evolved stars and could be formed within 
the interstellar medium (Hcnning, 199E). It is the moti- 
vation of this work to shed further light on the evolution 
of sma ll dust grains within the diffuse interstellar medium 
in rela ion to gas physical conditions. To do this we have 
combined ISOCAM observations of the Ursa Major cirrus 
cloud with 21 cm, CO and IRAS data. The Ursa Major 
cirrus has a moderate opacity to stellar light (Ay < 0.5) 
and an atomic and molecular sections traced by Hi and 
CO emission respectively. The comparison of these data 
sets allows us to measure the mid-IR emissivity per hydro- 
gen which, in a cloud transparent to stellar light, is related 
to the abundance and the optical properties of small dust 
particles independently of any modelling of the penetra- 
tion of the dust heating radiation. 

The observations used in this investigation are pre- 
sented in § ||. The structure and the kinematics of the 
Ursa Major cirrus is described in § [| In § [| we esti- 
mate the emissivity per H atom of small dust grains in 
the atomic part. Section || is devoted to the estimate of 
the abundance of small dust grains in an Hi filament with 
a particular kinematic and § [| addresses the same question 
for the H1-H2 interface. The relation between the spatial 
variations of the dust emission and the gas physical condi- 
tions are discussed in § [z| The main results of this analysis 
are summarized in § ^. 

2. Observations and data processing 

We present ISOCAM observations of th e Ursa Majo r cir- 
rus, a piece of the North Celestial Loop ( Heiles, 1989 ) (see 
Fig. ||). This cirrus is composed of two connected parts, 
one predominantly atomic, detected only in the Hi line 



Fig. 2. Field A : LW2 image (left), LW3 image (middle) 
and 21 cm integrated emission map (right). 



and one mostly molecular, detected in CO ( |de Vries et al 



1987; Pound & Goodman, 1997). For the sake of simplic 



ity we will name these two sections of the cirrus atomic 
and molecular in the following. The atomic section of the 
cirrus is representative of the so called Hi clouds with a 
density n ~ 100 cm -3 and a temperature T ~ 200 K 



(Joncas et al., 1992). It is illuminated by the mean solar 
neighbourhood interstellar radiation field (ISRF) and it 
is optically thin to this radiation; the Hi col umn densities 
range from 2 to 6 x 10 20 cm~ 2 . According to |de Vries et al 
(1987 ), the distance of the Ursa Major cirrus is ~100 pc. 



2.1. ISOCAM observations 

Three small fields of 0.05 square degree have been mapped 
with ISOCAM (see Table Q for the observation log). The 
location of the three fields is indicated in Fig. These 
fields have been selected according to the following cri- 
teria. First, field A crosses an Hi filament with a large 
transverse velocity gradient suggesting a rotation along 
the long axis (Miville-Deschenes, 1999). It provides the 
opportunity of studying the impact of this peculiar gas ve- 
locity distribution upon the spatial and size distributions 
of small dust particles. Secondly, fields B and C encom- 
pass an atomic-molecular interface, allowing the analysis 
of dust properties variations across this interface. 

The ISOCAM observations were made in the LW2 (5 
to 8.5 itm) and LW3 (12 to 18 /im) broad band filters. For 
diffuse clouds the emission in these filters is dominated 
by the zodiacal emission; the contrast of the cirrus is as 
small as a few 0.1 per cent in the LW3 filter. Furthermore, 
ISOCAM is affected by instrumental effects that can be 
one order of magnitude higher than the contrast we are 
looking for. To correct for these instrumental effects we 
have developed a data processing method based on the 
fact that a given point on the sky has been observed sev- 
eral times. T his method is described in Miville-Deschenes 
]pt al. (2000 ). Without this particular data processing, the 
ISOCAM observations presented in this paper would have 
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nela A 


nela d 


nela G 


Filters 


LW2 and LW3 


LW2 and LW3 


LW2 and LW3 


Right Ascension (J2000) 


9h 45m 24.8s 


9h 35m 7.4s 


9h 30m 45.6s 


Declination (J2000) 


70° 22' 54.8" 


70° 23' 15.1" 


70° 26' 48.4" 


Total number of readouts 


490 


487 


490 


Number of sky positions 


2 x 32 


2 x 32 


2 x 32 


Number of readouts per sky position 


5 


5 


5 


Step size between subsequent positions (pixels) 


8 


8 


8 


Integration time (seconds) 


5.04 


5.04 


5.04 


Lens (pixel field of view) 


6" 


6" 


6" 



Table 1. Description of the ISOCAM observations. The LW2 and LW3 broad filters correspond respectively to the 
wavelength ranges [5-8.5] and [12-18] /im. The size of the three fields is ~ 7' x 26'. 




9 h 3fl m 97 36 35 34 33 9 h 33^ 37 38 3$ 34 33 



Fig. 3. ISOCAM (LW2 and LW3), Hi and CO observations of field B. 




* h 33 m 32 31 30 29 2S B^aS™ 32 31 30 29 26 

Fig. 4. ISOCAM (LW2 and LW3), Hi and CO observations of field C. 
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Fig. 5. Column density profiles as a function of the dec- 
lination in field A for the Hi column density (black), LW2 
(blue), LW3 (yellow) and IRAS 100 /im (red). An offset 
has been subtracted to each profile (taken at the south- 
ern edge of the field). All the infrared emission profiles 
were normalized to the Hi column density by constant fac- 
tors: N LW2 = (/LW2 " 3.61)/6.5 x 10~ 22 , N hW3 = (/lws - 
19.12)/9.0 x 10~ 22 and iV 100 = (I wo - 2.7) /I A x lO" 20 



been useless. The LW2 and LW3 images of the three fields 
are shown in Figures |[ || and [I| The mid-infrared peak- 
to-peak fluctuations in these images are < 0.4 MJy sr _1 
and the rms of the high frequency noise is ~ 0.02 MJy 
sr _1 . 



2.2. Hi and CO observations 

We have combined the mid-infrared observations with CO 
and Hi observations to conduct our analysis. The Hi (21 
cm) observations were obtained at the Dominion Radio 
Astrophysical Observatory of Pcnticton (Canada). Part 
of these data were published by Joncas et al. (1992] ) and 
the whole set will be the subject of a forthcoming paper 



( [Miville-Deschenes, 1999| |Mivillc-Deschenes et al., 2001| ). 
The spatial resolution of these 21 cm observations is 1' 
and the spectral resolution is 0.66 km s _1 . 

In order to compare the dust emission with the whole 
gas content, we use 12 CO observations of the Ursa Major 
cirrus in the J=(l-0) transition at 115.271203 GHz that 
were obtained by Pound fc Goodman (1997| ). These au- 
thors kindly let us use their data obtained with a Gaussian 
beam of 100" and a velocity resolution of 0.27 km s _1 . The 
images of the integrated Hi and CO emission of the three 
ISOCAM fields are also shown in Figures |2[ || and £| 



3. General structure and kinematics 

3.1. Field A - an Hi filament 
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Fig. 6. Left: map of the 21 cm integrated flux in field 
A. Center: position (520oo)-velocity representation of the 
21 cm data. Right: three typical spectra with the result 
of the Gaussian decomposition. The positions of the three 
spectra are indicated with red points and dotted lines. 



3.1.1. Comparison of dust and gas tracers 

The field A observed by ISOCAM is a small cut across 
a filament (length ~ 5 pc and width ^0.5 pc), observed 
in all IRAS bands (see the 100 /im emission in Fig. [l]) 
and at 21 cm (Joncas et al., 1992). As seen in Fig. ||, the 
filament is also well detected in the mid-infrared. The CO 
observations of Pound fc Goodman (1997| ) do not include 
field A, but CO has been searched unsuccessfully at the 



0.1 K level within this area by de Vries et al. (1987) and 
G. Lagache (private communication). Therefore we will 
consider that field as predominantly atomic and use the 21 
cm data as a tracer of the spatial structure and kinematic 
of the gas. 

On a large scale, the structure in the LW2 and LW3 
broad band filters is very similar to the Hi emission. But a 
careful comparison of the spatial distribution of the mid- 
infrared and Hi emissions reveals significant discrepancies. 
In Fig. H we plot the total Hi column density, the LW2, 
LW3 and IRAS 100 ^im brightness, averaged over o^ooo? 
as a function of the declination 62000 in field A. In these 
profiles, the emission in the LW2 and LW3 filters and the 
100 /im emission have been normalized to the Hi column 
density with a constant ratio. These profiles show that the 
Hi emission is very well correlated to the 100 /im emission, 
but the peak of the emission in the two ISOCAM filters 
does not coincide with that of the FIR and Hi; the mid- 
IR emission is shifted compared to the Hi and 100 fim 
emissions. 

3.1.2. Kinematics of the gas 

The velocity structure of the Hi emission within field A is 
described in Fig. [©J. The central panel presents a velocity- 
position diagram of the 21 cm emission. Three represen- 
tative Hi spectra of this field are also shown in the right 
panel. On most of the field, the 21 cm spectra are charac- 
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Fig. 7. Upper panel: Mean integrated 21 cm emission as 
a function of the declination for the two Hi components 
identified by the Gaussian decomposition. Lower panel: 
mean velocity as a function of the declination for the two 
components. The continuous line depicts the filament and 
the dotted line depicts the more diffuse component. 



terized by two components. The two components clearly 
overlap in velocity and it is impossible to separate them 
only by selecting velocity ranges. To estimate the column 
density of each component across the field, we decompose 
the spectra by a sum of two Gaussians. From this de- 
composition, we derive the integrated intensity, the ve- 
locity centroid and the velocity dispersion of each com- 
ponent on the line of sight. The two components corre- 
spond to spatially coherent structures which extend be- 



yond the ISOCAM field ( [Mivillc-Dcschcnes et al., 2001| ). 
The brightest structure seen in the position-velocity dia- 
gram (Fig. H) at a velocity of ~ 8 km s~ 1 corresponds to 
the filament observed in the IRAS 100 fim map (Fig. |l|). 
The second component at a central velocity of ~ 4 km 
s~ 1 is more diffuse; it is also observed beyond field A. 

In Fig. |?] we show the mean integrated intensity and 
the velocity centroid of these two components across field 
A. The component associated with the filament (solid line 
in Fig. fj]) contributes ~ 40% of the total integrated emis- 
sion of the field. It is characterized by two velocity gra- 
dients perpendicular to the filament; first there is a 2 km 
s^ 1 increase of the velocity from £2000 = 70.15° to 70.35° 
and then a decrease of 4 km s _1 from ^2000 = 70.35° 
to 70.54°. At a distance of 100 pc, these variations cor- 
respond respectively to gradients of —5 and -10 km s _1 
pc -1 . The velocity dispersion of this component measured 
on individual spectra is ~ 2 — 3 km s . The more diffuse 
component (dashed line in Fig. 0), which contributes ~ 
60% of the gas emission, is characterized by a north-south 
integrated emission gradient and an almost constant ve- 
locity centroid (~ 4 km s _1 ). Its velocity dispersion is 
significantly larger than that of the filament, with values 
between ~ 6 and 10 km s _1 . 




143.5 



143,0 



142.5 



142.0 
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Fig. 8. Mean intensity in field C as a function of «20oo- 
From top to bottom: Hi(l'), CO (1.7'), IRAS 100 fira(4% 
IRAS 60 /im(3'), LW3 (6") and LW2 (6"). The numbers 
in parenthesis are the resolutions of the observations. 

3.2. Field B and C - an Hi -H2 interface 

The two other fields observed with ISOCAM cross an 
atomic/molecular gas interface as indicated by CO ob- 
servations. In field B (Fig. |J) there is a slight increase of 
the CO emission (up to Tb ~ 3 K km s _1 ) at the western 
edge of the field. The Hi emission in this field is very uni- 
form; the contrast in the 21 cm line integrated emission is 
-10%. 

In field C, Hi and CO observations are available for 
the whole region observed by ISOCAM. There is a slight 
north-south Hi emission gradient but again the contrast is 
very low (< 10%). Most of the column density variations 
in field C is associated with CO emission. 

In Fig. H we show cuts across field C in Hi, CO, IRAS 
100 fim, IRAS 60 /xm, LW2 and LW3. What is striking at 
first is the flatness of the Hi emission. There is a sharp 
rise of the CO emission near 02000 = 142.6°, not associ- 
ated with any change in the Hi emission. In contrast, the 
dust emission in all 4 bands increases in the molecular re- 
gion, but with less sharpness unrelated to the angular res- 
olution. The infrared dust emission systematically peaks 
on the western side of the CO filament. Furthermore, we 
also notice that the peak position of the dust emission at 



6 



Miville-Deschenes et al.: ISOCAM observations of the Ursa Major cirrus 



in 



14 ; 

12 ~- 

10 r 

8 

6 r 

4 r 
3 



'St. 



4.0 



a.e 



E 



2 r 



1 r 



L 




5 

v (km/s) 



10 



15 



Fig. 9. Top: 21 cm spectrum at the position of t 
CO J=(l-0) peak intensity in field C (02000 = Qh31m, 
^2000 = 70°23'). Bottom: CO spectrum at the same posi- 
tion. The result of a Gaussian decomposition of the spec- 
tra is also shown (dash lines). Two Hi components are 
identified (FWHM of 6.16 and 4.28 km s" 1 ). Only the 
largest of the two Hi components is seen in CO (FWHM 
of 1.50 km s" 1 ). 
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Fig. 10. LW2 (top) and LW3 (bottom) brightness as a 
function of the Hi column density (defined in the optically 
thin approximation). The linear regression over-plotted 
are: I hW2 (MJy/sr) = 6.07 x 10~ 22 N m (cm" 2 ) + 3.47 and 
hwa (MJy/sr) = 8.41 x 10~ 22 Af m (cm" 2 ) + 18.92 



shorter wavelength seems to be more shifted with respect 
to the CO emission than at longer wavelength. 

In velocity space, the Hi and CO data are very simple 
to describe. CO and 21 cm spectra taken at the peak CO 
intensity in field C are shown in Fig. |j. These spectra are 
representative of both fields B and C . The CO emission is 
characterized by a single symmetric component centered 
at ~ 1.3 km s _1 with a velocity width at half power of 
~ 0.7 km s _1 . The Hi emission is much broader; the spec- 
trum is clearly not symmetric with its peak intensity near 
~ 2 km s _1 . It can be fitted by two Gaussian components, 
one at v ~ 2.7 km s _1 with a width of ~ 2.7 km s _1 and 

a second one at v ~ 8.5 km s _1 with a width of ~ 2.0 km 

r,-i 



4. Emissivity of the small dust particles 

4.1. Correlation between the mid-infrared emission and 
the Hi integrated emission 

A main goal of this work is to study the relation between 
the dust and gas emission. If the dust-to-gas ratio and the 
heating of the dust are uniform and constant, the mid- 
infrared emission I m i r (a,S) may be described by: 



•(a,i) = ex N H (a,8) + hack 



(1) 



where e is the dust emissivity, Nh the total column den- 
sity of hydrogen (Hi and H2) and hack is the background 
level (assumed constant over the observed field). The total 
hydrogen column density is deduced from the in tegrated 
emiss ion at 21 cm (W m ): N m = 1.823 x 10 18 Wm flSpitzcr] 
1978| ) when the atomic gas is optically thin, which wc con- 



sider a reasonable assumption for this cloud. We did not 



look for a correlation between Njj and I m i r in fields B 
and C because of the very low contrast of the mid-infrared 
emission and because of the gre at uncertainty conce rning 
the CO to H2 conversion factor ( de Vries et al. (1987 ) give 
N H2 = 0.5 ± 0.3 x 10 20 Wco cm" 2 K" 1 (km/s)- 1 for the 
Ursa Major cirrus). 

The LW2 vs Nm and LW3 vs Nm plots in field A are 
shown in Fig. [l^. The mid-infrared emissions rise globally 
with the Hi column density; these relations are approxi- 
mately linear (see Fig. |l^). Nevertheless, the rms scatter 
in the LW2-iV Hl and LW3-7V H i relations (0.05 MJy/sr for 
LW2 and 0.07 MJy/sr for LW3) is significantly larger than 
the noise level (0.02 MJy/sr). 

4.2. Background emission in the ISOCAM images 

The scatter around the LW2 vs iVm and LW3 vs A^hi re- 
lations indicates that our description of the mid-infrared 
emission with Eq. [l] is probably too simple. This scatter 
may originate from a spatial variation of the emissivity or 
of the background emission. Here the background level is 
mainly due to the zodiacal light, which is uniform at the 
scale of our small fields. Therefore we conclude that the 
scatter of the LW2 vs Nm and LW3 vs A^hi relations in 
field A is due to spatial variations of the emissivity e. 

In order to quantify the mid-infrared emissivity for 
each position (a, 5), the uniform background emission 
level of the ISOCAM images must be estimated. This is 
also essential to compare the emission between the LW2 
and LW3 filters. We can estimate the background level 
by inspecting the relation between the mid-infrared fluxes 
and the total column density of hydrogen (see Fig [To]) . If 
the mid-infrared flux is, at first order, proportional to the 
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Fig. 11. Mean LW2 flux (at the resolution of the Hi obser- 
vation) as a function of the declination (continuous line). 
Left: the dotted line is the mid- infrared profile reproduced 
from the Hi data with an emissivity of e„ = 4.5 x 10 -22 
MJy/sr cm 2 (which corresponds to ^nve v — 3.8 x 10~ 31 
W/H) Right: the dotted line is a reconstruction of the mid- 
infrared profile with an emissivity of 8.2 x 10 -22 MJy/sr 
cm 2 (4.6 x 10~ 31 W/H) for the filament and of 1.6 x 1(T 22 
MJy/sr cm 2 (0.9 x 10~ 31 W/H) for the diffuse component. 
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Fig. 12. Mean LW3 flux (at the resolution of the Hi obser- 
vation) as a function of the declination (continuous line). 
Left: the dotted line is the mid-infrared profile reproduced 
from the Hi data with an emissivity of 7.1 x 1CP 22 MJy/sr 
cm 2 (2.3 x 1(T 31 W/H). Right: the dotted line is a recon- 
struction of the mid-infrared profile with an emissivity of 
1.1 x 1CT 21 MJy/sr cm 2 (2.5 x 1CT 31 W/H) for the fila- 
ment and of 3.8 x 1(T 22 MJy/sr cm 2 (1.0 x 10~ 31 W/H) 
for the diffuse component. 



total column density of hydrogen, the background level 
should be the mid-infrared flux for zero column density of 
hydrogen. If we apply this to the data the uncertainty on 
the background will be large because the dynamical range 
of the ISOCAM images is small (< 0.6 MJy/sr). Further, 
in fields B and C, CO emission cannot be considered as a 
quantitative tracer of the H2 column density. In this con- 
text, we chose to work in relative brightness for all tracers. 
This means that, for each field and for gas and dust maps, 
we subtract an offset taken in a common region where the 
mid-IR emission is the lowest. 

5. Variation of the mid-infrared emissivity with 
gas velocity 

In this section we show that the spatial shift between the 
dust and gas tracers in field A, and therefore the dispersion 
in the LW2-./Vhi and LW3-iVni relations, can be accounted 
for by a difference in the small dust particles emissivity 
per H atom between the two Hi kinematic components 
(see § |3.1.2|). 



As seen in Fig. [?], the two Hi components identified ear- 
lier with the Gaussian decomposition have very different 
column density profiles across field A. If we consider the 
morphology of both components, the mid-infrared emis- 
sion is more similar to the filament than to the diffuse 
Hi component. By giving different emissivity to these two 
Hi components we were able to reproduce rather well the 
mid- infrared profiles (see Fig [n] and |l^). The LW2 and 
LW3 emissivities that best reproduce the mid-infrared ob- 
servations are detailed in Table ||. 



The comparison of the mean emissivity of typical 
Galactic cirrus of the solar neighborhood obtained at 12 
^m by |Boulangcr fc Pcrault (1988] ) with the LW2 and LW3 
emissivities (Airve^ in W/H) deduced here, shows that the 
diffuse component of field A has an emissivity very sim- 
ilar to the average cirrus value. On the other hand, the 
filamentary component seems to have an emissivity 5 ± 1 
(LW2) and 2.5 ± 0.5 (LW3) times higher than the average 
value. From the numbers of Table || we also conclude that 
the ratio of the LW2 and LW3 emissivity is ~ 0.9 in the 
diffuse gas and ~ 1.8 in the filament. Because of the good 
spatial correlation between the 100 fim and Hi column 
density (see Fig. ||), the large grains responsible for the 
100 /im emission must have the same emissivity in the the 
filamentary and diffuse components. 

6. Abundance variations of small dust grains in an 
H1-H2 interface 

Fields B and C cross an interface between a mainly atomic 
region and another one where significant CO emission is 
detected. The emission in field B is very uniform: it is thus 
difficult to correlate variations in the dust and gas emis- 
sions. Therefore, our whole analysis of the H1-H2 transi- 
tion will be based on the observations of field C . 

To estimate the dust mid-infrared emissivity, we usu- 
ally rely on an estimate of the gas column density. In field 
C, the structure of the gas is dominated by the molecular 
component. Because the determination of the CO to H2 
ratio is poorly constrained and because, in such diffuse 
clouds, it is known that CO is not such a good tracer of 
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n q -i_ n q 
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filament 
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2.3 ± 0.1 
2.5 ± 0.3 
1.0 ±0.3 


(K km/s) 


c 


6.5 (LW2) 


Hi 

CO 


1.9 ±0.2 
0.48 ± 0.05 


S 




15.0 (LW3) 


Hi 

CO 


1.3 ±0.1 
0.61 ±0.06 




high latitude 


12.0 




1.1 ±0.3 


o 




100.0 




3.2 ±0.1 


o 



Table 2. Emissivity in LW2 and LW3 in fields A and 
C. The 12 and 100 /im emissivity of typical h igh lati- 
tude Galactic cirrus ( Boulanger fc Pcrault, 1988 ) are also 
given as a comparison. The uncertainties on the emissiv- 
ity for the diffuse and filamentary components of field A 
reflect the uncertainties we have on the determination of 
the background emission. The uncertainty on the average 
emissivity of field A has been computed using the boot- 
strap method flEfron fc Tibishirani, 1986; ) on the LW2-N ffi 
and LW3-Nhi relations of Fig. |lGj. The uncertainty on the 
emissivities of field C reflects the variations of the emis- 
sivity values over the region where they were measured. 




143 



142,8 



142.6 142.4 
g^o (degree) 



142.2 



142.0 



the molecular column density (Reach et al., 1994), we pre- 
fer to compare the mid-infrared emission with the IRAS 
100 /im band emission to estimate the abundance of the 
smallest dust particles, relative to that of bigger grains 
responsible for the emission at 100 /im. 

In Fig. |l3| (middle panel) we show the ratio 
(yl v )x/(ul u ) [iaaiim) for LW2 (A=6.5 /mi), LW3 (A=15 
/mi) and IRAS 60 /mi, across the H1-H2 transition. This 
way of presenting the data emphasizes the variation of 
dust emissivity at wavelength A compared to the emis- 
sivity at 100 /im. If the grain properties and the heating 
are homogeneous, this ratio reflects changes in the rela- 
tive abundance of grains. For all three wavelengths a de- 
crease in the (y 1 v ) \ / \v 1 ~ v ) ( 100 ratio is observed, corre- 
lated with the increase of the CO emission. Furthermore, 
for the three bands, the lowest (vIu)\/(vIv)(ioonm) is ob- 
served near the CO peak emission. 

By making the assumption that the 100 /im emissiv- 
ity is that measured by Boulanger & Perault (1988) for 
the high latitude clouds (see Table |2|), which is in ac- 
cordance with the 100 /im emissivity estimated in that 
cloud by Joncas et al. (19*9^ ), we were able to estimate 
the emissivity in the LW2 and LW3 bands from the 
(yl v )\j {vlv)(i00fj,m) ratio. In Table ||we list an estimate of 
the mean LW2 and LW3 emissivity values for the atomic 
and molecular regions of the H1-H2 transition. We esti- 
mate that the ratio of the emissivity in the Hi region com- 
pared to the CO region is 4.0 ± 0.8 and 2.1 ± 0.4 in LW2 
and LW3 respectively. 



Fig. 13. Top: Mean CO emission of the western 
part of field C, as a function of 0:2000 • Middle: 
Ratio (7 hu2 /6.5)/(/ioo/100), (7 to3 /15)/(/ioo/100) and 
(/6o/60)/(/ioo/100) as a function of 02000- Bottom: Ratio 
of the LW2 over LW3 emission as a function of 0-2000- 

The ratio of the LW2 over LW3 emission for field C 
is also shown in Fig |l^. In the purely atomic part of the 
fields, the ratio is ~ 0.8 and it drops to ~ 0.3, in a few 
O.lpc, as we enter the molecular region. 

7. Discussion 

The comparison of the mid-infrared, Hi, CO and IRAS 
observations has allowed us to highlight important spatial 
variations of the small dust particles emissivity in the Ursa 
Major cirrus. 

7.1. Spatial variation of the mid-infrared emissivity 

The Ursa Major cirrus is characterized by a low gas col- 
umn density (6.0 x 10 20 cm -2 - corresponding to Ay ~ 0.3) 
and the absence of star formation. In this context, the 
spatial variations in the mid-infrared emissivity can re- 
sult from either (1) the presence of molecular gas or very 
cold Hi gas or (2) an attenuation of the UV radiation 
field on scales smaller than the beam size or (3) varia- 
tions in the small grains abundance or emission proper- 
ties. The good correlation observed between the Hi and 
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100/im emissions in Fig. ^| puts constraints on hypothesis 
(1) and (2). Dust associated with molecular gas should 
produce an excess of IR emission with respect to the IR- 
Hi correlation. However, this excess could be reduced by 
an attenuation of the UV radiation field in the same re- 
gions since molecules are expected to be present in the 
most shielded parts of the cloud. But both effects would 
modify the mid-infrared and 100 /im emission in the same 
way. Therefore, variations of the dust size distribution or 
of emission properties (hypothesis (3)) are more directly 
traced by the mid-infrared to 100 /im emission ratio than 
the mid-infrared emissivity per H atom. However it is only 
with the Hi emission that we have access to the velocity 
information. 

In field A, the emissivities in the LW2 and LW3 filters 
are enhanced by a factor 5±1 and 2.5 ±0.5 in the filament 
with respect to the more diffuse atomic gas. Since such a 
change is ruled out for large grain emission at 100 /im, we 
interpret the LW2 and LW3 emissivity variations as an 
evidence for variations in the small grains characteristics 
(abundance/emission properties). Elements are still lack- 
ing to identify whether the emissivity variations oberved 
here are due to variations in the abundance of the mid- 
infrared emitters or to their emission properties. We re- 
call that this filament has a very peculiar kinematic: a 
large velocity gradient across the filament which could be 
a signature of a fast rotation around its long axis. We dis- 
cuss below how large velocity gradients might enhance the 
abundance of small dust grains. 

Based on the mid-infrared to 100 /im ratio, in field 
C, the abundance of the particles emitting in the LW2 
and LW3 filters are found to be respectively 4.0 ± 0.8 and 
2.1 ± 0.4 times lower in the molecular region than in the 
atomic gas. A similar effect is also observed with the IRAS 
60 /im emission. 

7.2. Turbulence and dust size distribution 

We discuss here a possible link between the dust size distri- 
bution and the kinematics of the gas (i.e. the importance 
of the turbulent motions). The filament in field A, which 



has an aspect ratio of the order of 10 ( Miville-Deschenes 



199E), is reminiscent of coherent structures of vorticity 



observed in turbulence and thought to be responsible for 



the phenomenon of intcrmittency (Vincent & Meneguzzi 
1994). Falgarone fc Pugct (1995| ) investigated the decou- 



pling of dust and gas in a turbulent flow as a function of 
grain size. They showed that the intermittency of turbu- 
lence enhances the decoupling of small dust grains from 
the gas motions because of the non-Gaussian distribution 
of the velocity gradients. Thus the relative velocities of 
small dust particles can be significantly increased above 
that reached in a Kolmogorov turbulence. Qualitatively, 
the collision outcome is coagulation or grain fragmentation 
below and above some energy threshold. In these regions, 
the relative velocities might become sufficiently large that 
small dust grains are produced by grain fragmentation in 



grain-grain collisions. The high abundance of small dust 
particles in the filament of field A could be a signature 
of such a fragmentation process. In the molecular gas, 
where the turbulent motions are smaller than in the Hi 
gas (Fig. ^), we may observe the reverse process: mutual 
coagulation of small grains could occur which would ex- 
plain the drop of the mid-infrared emissivity. 

The outcome of a grain-grain collision (coagulation, 
bouncing or fragmentation) depends on the kinetic energy 
of the colliding grains, on their respective masses and on 
their material. Based on a numerical simulation of aggre- 
gate growth, Dominik & Tielens (1997) have showed that 
two ~ 0.1 /xm aggregates stick if their relative velocity is 
8v < 0.8 km s _1 . But if it is higher than ~ 1 km s -1 , the 
dust aggregates start to fragment. 

Following Draine (1985), the typical relative velocity 
of two grains of size eji in a Kolmogorov turbulent flow 
characterized by a velocity dispersion v max at a scale l max 
is given by the following equation: 



Vrel 




l/ 2 / 2tt \ 



1/4 



(2) 



where n and T are the gas density and temperature and 
p gr is the grain mass density. If we consider standard 
values for an Hi cloud like Ursa Major (n = 100 cm -3 , 
T = 200 K) and Kolmogorov turbulence with l max = 10 



pc and Vn 



5 km s , the typical relative velocity for 



big grains (a± = 0.1 /im and p gr — 3 g cm 3 ) would be 
v re i ~ 0.6 km s _1 , which is close to the fragmentation 



threshold established by Dominik fc Tielens (1997| ). But 
in field A, we observe velocity gradients of 5 and 10 km 
s _1 pc -1 , which are larger than the gradients characteris- 
tic of Kolmogorov turbulence at the scale of the filament 
(0.3 pc). As shown by Falgarone fc Puget (1995| ), these 



local large velocity gradients may increase the relative ve- 
locity between small grains. Considering the strength of 
the velocity gradient observed here, it is very likely that 
the relative velocity between large grains exceeds the frag- 
mentation threshold. 

In field C, the width of the CO line is ~ 1.5 km s _1 
(see Fig. ^). If we consider that this width is mostly due 
to turbulence and that the thickness of the Ursa Major 
cirrus is ~ 2 pc, the typical velocity dispersion of 0.1 /im 
grains for the same temperature and density as in field 
A is ~ 0.2 km s _1 . This is lower than the coagulation 



threshold determined by Dominik & Tielens (1997) and 



it could be even lower since the density in the molecular 
region is higher than the mean cirrus value. 

This indicates that fragmentation and coagulation by 
grain-grain collisions in a turbulent flow may be at the 
origin of the abundance variations of small dust grains ob- 
served in the Ursa Major cirrus. To estimate the timescales 
over which the dust size distribution is significantly af- 
fected by this process, one should take into account the 
collisions between grains of all sizes and, in the case of 
fragmentation, make some hypothesis on the size distri- 
bution of the fragments. This is beyond the scope of the 
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present paper. Drainc (1985) has estimated the deple- 
tion time scale of small grains in interstellar clouds to be 
~ 10 6 years by considering only sticking on large grains 
in Kolmogorov turbulence. Mutual coagulation of small 
grains actually occurs faster if we take intermittency into 
account (Falgarone & Puget, 1995). 



7.3. Variations of the I U (LW2)/I U (LW3) ratio 

In addition to the spatial variations of the mid-infrared 
emissivity we have also observed spatial variations of the 
J„(LW2)//„(LW3) ratio. For example, it is twice larger 
in the Hi filament of field A than in the diffuse gas that 
surrounds it. It is also observed to be twice larger in the 
atomic side of the H1-H2 transition than in the molecular 
region. Similar variations of the J„(LW2)/I„(LW3) ratio 
are observed across photo-dissociation regions from the 
photo-dissociated gas to the molecular gas ( |Abergel et al 
200C| ). 



Usi ng broad-band filter s and CVF observations of 
L1630, |Abergel et al. (2000|) showed that the emission 111 
the LW2 filter is correlated to the intensity of the 7.7 /im 
aromatic band, and that the emission in the LW3 filter is 
correlated to the continuum level at 15 /im. These CVF 
data show that variations in the relative amplitude of the 
aromatic bands are minor and do not significantly affect 
the J l/ (LW2)/J ly (LW3). It is unclear whether the aromatic 
bands and continuum at 15 /im come from the same dust 
particles. Within the present understanding of the nature 
and emission properties of small dust grains it is impossi- 
ble to determine the origin of the i !/ (LW2)//„(LW3) ratio 
variations. 

7.4. Coagulation of small grains 

Another example of the decrease of the small dust grain 



resolution and brightness sensitivity of these images im- 
prove by one order of magnitude those provided by the 12 
and 25 /im IRAS images. We have compared the ISOCAM 
images with Hi, CO and 100/im observations to mea- 
sure the mid-IR emissivity (mid-IR emission from small 
dust particles per hydrogen atom) for distinct gas veloc- 
ity components and across an interface between atomic 
and molecular gas. We find significant variations in the 
mid-IR emissivity which we interpret as changes in the 
abundance of small grains within the cloud. The abun- 
dance of the small dust particles is larger by a factor 5 
within a filament with a strong transverse velocity gradi- 
ent, than in the surrounding diffuse Hi. We propose that 
this abundance enhancement results from the formation 
of small dust grains through the fragmentation of larger 
grains in energetic grain-grain collisions induced by turbu- 
lence. The abundance of small dust particles with respect 
to the larger grains is observed to be reduced (by a fac- 
tor 4) in CO emitting molecular gas with respect to its 
value in diffuse Hi. In this gas the turbulent motions have 
a smaller amplitude than in the Hi gas. Within the pro- 
posed interpretation, we speculate that in such an envi- 
ronment grain/grain collisions are less energetic and they 
therefore lead to coagulation and not fragmentation. 

Spectroscopic observations indicate that the 
LW2/LW3 color ratio measures the contrast of the 
aromatic bands with respect to the underlying continuum 
emission. A drop in the / ! ,(LW2)//„(LW3) ratio by a 
factor — ■ 2 is observed from the Hi filament to the diffuse 
Hi gas and from the Hi to the H2 gas across the Hi 
-H2 interface. The LW2/LW3 color ratio is observed to 
decrease with the mid-IR emissivity. If the band con- 
trast decreases for increasing particles sizes, a plausible 
assumption in terms of solid state physics, the mid-IR 
emissivity and color variations could both be related to a 
change in the small-sizes end of the dust size distribution. 



abundance in an H1-H2 interface is given by Bernard et al. 
1999| ) who characterize dust properties, in a section of Acknowledgements 



the Polaris cirrus with moderate opacity (Ay ~ 1), by 
combining IRAS and far-IR/sub-mm observations. In this 
cloud, there is no mid-IR emission, the 60/100 /im emis- 
sion ratio is low and the large grain temperature, 13 K, 
is lower than the mean cirrus value 17.5 K. These dust 
properties are characteristic of dense gas within molecu- 



lar clouds. Bernard et al. (1999) explain these observations 
by the formation of large fractal grains through the coag- 
ulation of small dust grains on large grains. The decrease 
of the small dust grains abundance in the H1-H2 interface 
of the Ursa Major cirrus could be an intermediate step 
in the dust evolution, between the diffuse Hi and denser 
environments. 



8. Conclusion 

We have presented ISOCAM images (in two filters LW2 
(5-8.5 ^m ) and LW3 (12-18 /im )) of a high Galactic lat- 
itude cirrus transparent to stellar photons. The angular 
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